SOFT INTERACTIONS 



o 



X 



RALPH ENGEL 

Bartol Research Institute, University of Delaware 
217 Sharp Lab, Newark DE 19716, USA 
E-mail: eng&lepton. bartol. udel. edu 

A brief introduction to the theory and phenomenology of soft interactions is given, 
focusing on total and elastic cross sections and multiparticle production. 



1 Introduction 
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\^ , The description of soft hadronic multiparticle production within QCD is one 

of the long-standing, unsolved problems. In recent years progress.^ has been 
. made in understanding soft QCD processes in lattice calculationstJ and also 

' the Jiodified leading-log approximation proved to be a remarkable powerful 

Q>^ . toolo. In this article we focus on another method which is the investigation 

CO ' of predictions implied by unitarity and analyticity of scattering amplitudes 

and general features of gauge field theories. We will discuss basic concepts 
and predictions of general nature rather than one or several models in detail, 
emphasizing a pedagogical presentation of the underlying ideas. The aim is 
to supply the reader with a basic theoretical framework to understand the 
more specific contributions to this conference. More detailed discuffikma of 
soft interactions and related topics can be found in other reviews 



2 Regge amplitude 



Cross sections and amplitudes at low energy are fairly well understood in terms 
' of production and decay of a single or a few resonances, as shown in Fig. 

With increasing collision energy more and more resonances (with increasing 
mass and decay width) contribute to the total cross section. This seems to 
render any approach, which is based on summing all resonance contributions 
to extrapolate to high energy, infeasible. However, under certain conditions 
this summation can be carried out and the Reggeon amplitude is obtained u. 

Because of angular momentum conservation it is convenient to expand the 
scattering amplitude A(s,t) in terms of partial waves a;(s) for fixed angular 
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(a) (b) 

Figure 1. Production and decay of A++ in the s-channel (a), and a interaction obtained 
by crossing with A++ in the t-channel (b). 



momentum I 

oo 

A{s,t) = 16n^{2l + l)ai{s)Pi{cose), (1) 

(=0 

with 9 being the scattering angle in cm. frame and P; the Legendre functions 
of the first kind^. For resonance exchange the amphtude is of the weU known 
Breit-Wigner form a; ^ l/(s — mf + imiTi). 

One of the important assumptions in Regge theory is that of maximum 
analyticity of the scattering amphtude, which holds in gauge field theories. 
The amplitude is an analytic function and has only physics-implied singulari- 
ties. Thus, if the amplitude is completely known for a given scattering process, 
then the amplitudes for all other channels, which are related to the former 
one by crossing, can be obtained by analytic continuation of this amplitude. 
An example for this is given in Fig. 0b. After crossing from the s-channel to 
the t-channel Eq. (|^) reads 

A{s,t)^16TTY^{2l + l)ai{t)Pi{zt), zt=cos0t = -^ + l . (2) 

The scale Sq is related to the masses of the particles (sq ~ IGeV^). The 
partial wave amplitude describes now the exchange of a particle with angular 
momentum I in the f-channel ai{t) ^ l/{t — mf + imiTi). 

Using Cauchy's theorem the sum (|^) is re- written as integral (Sommerfeld- 
Watson transformation) 

Ms,t)=J2^[ ^U2^ + l)f ^t7r' ) ai{t) P,{~zt), r = ±l, 
r~±i 2* ■'Ci \ sm(7r0 J 

(3) 

"For sake of clarity we omit here all spin-related complications. 
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where the signature r = ±1 separates the sum into even and odd integer 
angular momenta. The sphtting of the partial sum into positive and negative 
signature contributions is needed to ensure convergence of the integral. The 
integration contour runs along both sides of the positive 5Re{Z} axis in the 
complex I plane (see Fig. ||a). 

In (||) the partial wave amplitude has to be extended to continuous values 
of I. This is done by employing the experimentally established, approximate 
relation between resonance mass and total angular momentum mf = al + ttiq 
with a and mo being constants specific to the group of resonances. Fig. ^ 
shows an example for meson resonances. From the structure of the partial 
wave amplitude 

1 ^ 1 1 ^ 1 

t'^ — mf t — itLq ~ al l — t/a + m'^/a I — a{t) 

follows that there is a singularity for I = a{t) with a{t) — (t— ttiq) / a. The func- 
tion a{t), called Regge trajectory, is in general a nonlinear complex function 
of t and related to the quantum numbers of the hadrons whose mass-angular 
momentum relation it describes. 




012345678 
(a) (b) spin I 



Figure 2. (a) Integration contours for Sommerfeld-Watson transform, and (b) mass-spin 
relation for for p and s mesons (Chew-Frautschi plot). 



Assuming that / = a(t) is the only singularity due to ai{t) we can dis- 
place the integration contour in Fig. ||a from Ci to C2 . The behaviour of the 
integrand leads to a vanishing contribution from the semi-circle at infinity. 
Neglecting the contribution from 3m{Z} = —1/2 all that is left is the con- 
tribution of the pole at I = a{t). Thus the final amplitude can be written 
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as 

1 I -i7va{t) 

sin(7ra(t)) ^ ' 

with I3(t) standing for the residue and other unimportant factors. 

Finally, the Regge amplitude is obtained by going to the high-energy limit 
s > So 

2s \ / s ^ 

■^-^ - \ s — >oo I ^ 



and reads 

A(s,t) = ^r;K(i))/3,(t)(^^y"^*^ . (7) 

The index k represents the sum over all contributing Regge trajectories and 
ri is the signature factor of the Regge trajectory k 

ri{au{t))^~—-^—r-. (8) 
sm(7rafc(t)) 

A single term in the sum (j^) describes the contribution of an infinite number 
of particles exchanged in the i-channel. Each of these groups of particles 
is represented in the amplitude by a quasi-particle, called reggeon, with the 
same quantum numbers. 



3 Cross section phenomenology 

3. 1 Regge- based fits 

Knowing the Regge trajectories from the measured masses of mesons and 
baryons, Eq. (0) can be used to predict the energy-dependence of elastic and 
total cross sections 

^ = ^\A{s,t)\\ a,ot^-^m{A{s,t = 0)). (9) 

dt iDTTS^ S 

The last relation employs the optical theorem and follows from unitarity^. 
Using furthermore the parametrization f3'^{t) = g^expjBoi} and considering 



''The masses of the scattering particles are neglected. 
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only the highest Regge trajectory, one gets with A — q;(0) — 1 

2A 



Cola 



CTtot = g 



■So 



9 

IGtt 

A 



So 



exp {Bolat} , f. 
la = Bo + 2a' (0) Ins 



3to A(s, t) 



t=o 



(10) 



There are several shortcomings of the amplitude ^ . The total cross sec- 
tion is predicted to have the energy dependence crtot ~ ■ The intercept of 
the highest known Regge trajectory is about a(0) « 0.5. The experimentally 
observed rise of the cross sections at high energy has lead Pomeranchuk to pro- 
pose the existence of a pole with a(0) ^ 1 which corresponds to the exchange 
of objects with vacuum quantum numbers. The corresponding quasi-particle 
is called the pomeron. Its trajectory should be related to glueballs instead of 
mesons or baryons. So far the existence of glueballs could not be confirmed 
experimentally. 

Another problem of Eq. is the violation of unitarity at high energy. 
This is most obvious by noting that the elastic cross section is predicted to 
finally exceed the total one or that the Froissart-Martin bound (T*°* < c In^ s 
is violated. 

Despite of these shortcomings, Regge pole motivated parametrizations 
are rather successful in pmdicting high-energy cross sections in the currently 
accessible energy range For example. Fig. ^ja shows the prediction by 
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Figure 3. Comparison of a simple Regge fit using-itwo poles with pp and pp data. The 
parameters for the total cross section fit are from t3. 

Donnachie and Landshoff together with pp and pp data. In Fig. we 
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show the energy dependence of the elastic slope parameter, Beiaj which is also 
in agreement with the simple Regge pole expectation of a logarithmic energy 
dependence (see Eq. (px|)). 

3.2 Multi-pomeron models 

As already mentioned, the pomeron exchange amplitude violates unitarity. 
Although this doesn't seem to be a serious problem at current collider energies, 
it indicates an inconsistency of this model. 

The simplest way of overcoming this problem without abandoning Regge 
theory altogether is the introduction of multiple pomeron exchanges in a single 
scattering process. Assuming that large momentum transfers are suppressed 
by the dynamics of the strong interaction, as it is the case in QCD, multi- 
pomeron exchange can be described on the basis of Gribov's Reggeon field 
theory (RFT) H. Then the total amplitude can be written as the sum of 
n-pomeron exchange amplitudes A^'^\s,t). For each n-pomeron graph one 
can define a theoretical "total" cross section applying the optical theorem to 
the corresponding n-pomeron amplitude 

oc 

(t(") = (-l)"+i-5m , (Ttot = ^(-1)"+V(") . (11) 

Here the alternating sign has been introduced by definition to keep all partial 
cross sections cr*^"^ positive. 

As a simplified model we consider only the first two graphs shown in Fig. ^, 
assuming ct'-"-' ^ < (T*-^-' with n > 2. Then, the total cross section reads 



(a) (b) (c) 

Figure 4. Hadron-hadron scattering via pomeron exchange: (a) one-pomeron, (b) two- 
pomeron, and (c) three-pomeron exchange graphs. 



Ctot = o'^^'^ — (7^'^\ where cr^^^ and cr'^^ are the cross sections of the one- and 
two-pomeron exchange graphs, respectively. The energy-dependence of the 
two-pomeron cross section is directly linked to that of cr'-^-' ~ and turns 
out to be cr*-^^ ~ s^^ . The two-pomeron cross section grows faster with energy 
than the one-pomeron cross section. This leads to a weaker energy-dependence 
of the total cross section than in the single-pomeron exchange model. 
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One particular realization of the multi-pomeron exchange is the eikonal 
model using the pomeron amplitude as eikonal function. In general, pre- 
dictions of such models are in good agreement with cross section data and 
are nearly indistinguishable from the single-pomeron exchange model in the 
energy range up to several TeV c.m.s. energy. 



3.3 Photon-induced hadronic interactions 

The generalization to photon-induced interactions is straight-forward in terms 
of the generalized vector dominance model (GVDM) E3o. The photon is con- 
sidered as a superposition of two types of states, a bare photon and hadronic 
fluctuations. With the probability of the order of the fine structure constant 
Qfem the photon can be found in a fluctuation which interacts like a hadron. 
Since the cross section of the bare photon is small it can be neglected in 
most applications. In the proton rest frame, the life time tf of the hadronic 
fluctuation with momentum k and mass M is much larger than the typical 
interaction time ti 

1 _ 2k 

^ VM2 + /fc2 _ ^_Q2 + p M2 + Q2 ' y ) 

where denotes the photon virtuality. The hadronic amplitudes for the 
scattering of longitudinally (i) and transversely (T) polarized photons are 
written as 

with g, my, and Tv being the four- momentum of the incoming photon and 
the mass and the decay width of the vector meson V , respectively. Ayf^^-^ is 
the amplitude for the scattering of the transversely polarized vector meson V 
off the hadron h. The parameter is introduced to allow for different cross 
sections of transversely and longitudinally polarized vector mesons. 

The GVDM is very successful in predicting_the total and and quasi-elastic 
vector meson photoproduction cross sections Ej. However, the s-dependence 
of the 7p and 77 cross sections changes considerably from real to virtual 
photons. This cannot be understood on the basis of GVDM alone. Many 
alternative models have been published. The most promising calculations 
combine within the color-dipole framework ideas of GVDM with perturbathie 
QCD calculations, obtaining a reasonable description of the measurements t3. 
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3-4 Impact parameter picture of high-energy scattering 

The impact parameter picture of high-energy scattering is very useful for 
discussing effects related to the wave character of the scattering particles. 
Its meaning becomes most transparent by considering again the partial wave 
expansion of the scattering amplitude 

oo 

A{s,t) = 167r^(2? + 1) ai{s) P;(cos6l). (15) 

;=o 

In the classical (geometrical) description, the ^th term on the r.h.s. of Eq. ( |l5|) 
with angular momentum / is assumed to correspond to the interaction of two 
particles with an impact parameter b satisfying I = \b\k. The impact parame- 
ter vector b is by definition perpendicular to the momentum k of the incoming 
particle. At high energy the sum can be approximated by the integral 

/•oo 

A{s,t)^167T dl (21 + 1) ai{s) Pi{cos0) , (16) 

"'0 



/ — ^oo 



using for a/(s) the analytic continuation in I. Taking the limit P;(cos( 
Jo [(2/ + 1) sin(6'/2)] and substituting 

Jo{z) = ^ f d^e'''°'V, b-q = b-q^=bq^cos<f , (17) 
271- Jo 

Eq. ( p^ ) becomes 

A{s,t) ^ As I d^b a{s,b)e'^^-\ with a{s,b) ^ al{s)\^^^g . (18) 



The impact parameter amplitude a(s, b) is the Fourier transform of the elastic 
scattering amplitude. Furthermore, in analogy to geometrical optics, the func- 
tion a(s, b) can be interpreted as the density function of sources of scattered 
waves producing interference patterns. The total and elastic cross sections 
are given by 

^Tcia(s) = 4 J d^b\ais,b)f, fTtot(s) =4 J d^b'imiaisM) . (19) 

For an amplitude with exponential t-dependence A{s, t) ~ 
istJtot exp{ii?ciat} the corresponding impact parameter amplitude reads 
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Figure 5. Partonic inter- 
pretation of the logarith- 
mic growth of hadrons with 
energy. 

gray" area) is given by the 
t-slope of the elastic cross section and, according to Eq. grows logarith- 
mically with s. This can also be understood in terms of partons in a hadron 
as shown in Fig. ^ Each parton splitting involves a transverse momentum 
transfer which corresponds to a small displacement in impact parameter. The 
total displacement of a parton follows from the number of steps N of its ran- 
dom walk, ^ Ab ■ N {Ap±)^^ Ins. This picture also readily explains 
why the elastic slope of hard processes has a much weaker energy dependence: 
the displacement in b is much smaller because of the larger Ap± . 

From the unitarity constraint |a(s,6)| < 1 one gets as high-energy limit 
ftot < cln^s with c being a constant. The special case of the black disk 
limit (e.g. maximum absorption) corresponds to a(s,6) = i/2 for 5 < i? 
and a(s, b) = otherwise. Then the inelastic cross section (cross section for 
absorption) is given by the geometrical size of the disk tJine = Tri?^ whereas 
the total cross section is twice the disk size, (Xtot = 27ri?^. This is a result of 
unitarity: absorption gives rise to elastic scattering a^ia = t^R^ = Ctot/2. 

4 Multiparticle production 

4-1 Unitarity and cuts 

To link particle production to elastic scattering amplitudes one can use the op- 
tical theorem. In terms of Feynman diagrams a particle propagator collapses 
to an external on-shell particle line if its imaginary part is taken 

3m ( 5 -] = tt5{p^ — m^) 

Therefore taking the discontinuity (i.e. the imaginary part) of a propagator 
is commonly referred to as "cutting" . This is also the reason why, up to kine- 
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^Pt) ^B) ~ 1 

R~log(s) 

The width of the amplitude in 6-space (i.e. the " 



(21) 



matical factors, the discontinuity of elastic amplitude is equal to the squared 
matrix element describing particle production, as sketched in Fig. 0. 



Figure 6. Graphical 
representation of 
optical theorem. 



One particular application is the discussion of the topologies of cut 
pomeron and reggeon amplitudes in the limit of large numbers of colors and 
flavours in QCD Ell. As depicted in Fig. Q the dominant final state contribu- 
tions are a two- and one-string configurations. 




(a) (b) 



Figure 7. Cuts and final state topologies of (a) pomeron and (b) reggeon exchange graphs. 



4-. 2 Abramovski-Gribov-Kancheli cutting rules 

In general, there are many possibilities of cutting diagrams involving multi- 
ple pomerons. The dominant contribution of the different cut configurations 
to the total discontinuity can—be calculated with the Abramovski-Gribov- 
Kancheh (AGK) cutting rules M: 

(i) The dominant contributions are given by cut configurations where the cut 
involves all intermediate particle states of a pomeron (a contribution due to 
a partially cut pomeron is sub- leading) . 

(ii) For an 7i-pomeron exchange graph, the contribution to the discontinuity 
with v cut pomeron propagators is B" 25m Here, A^"' denotes the 
n-pomeron amplitude and v is restricted to < < ^- The combinatorial 
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factors (AGK weights) are 



, ^ ' v\[n — v)\ - , , 

B:={ (22) 

1 - 2"-i : 1/ = . 

(iii) The coefficients -B" satisfy X]"=o-^" = 1' which means that all leading 
contributions to the total discontinuity are included. 

For example, there are three leading cuts of the two-pomeron graph 
(Fig. ^p): the diffractive cut with the weight -1 (Fig. ||a), the one-pomeron 
cut with the weight 4 (Fig. ^), and the two-pomeron cut with the weight -2 
(Fig. ||c). The cross sections for the different final states of the two-pomeron 



(a) (b) (c) 

Figure 8. Breakdown of the total discontinuity of the two-pomeron exchange graph ac- 
cording to the AGK cutting rules: (a) the diffractive cut describing low-mass diffraction, 
(b) the one-pomcron cut, and (c) the two-pomeron cut. 

exchange graph are (a) diffractive cut: (Tdm ~ ^^"^K (b) onc-pomcron cut: 
ui = — 4cr'^\ and (c) two-pomeron cut: 02 = 2a'^'^\ see Fig. ^. Note that the 
diffractive cut of the two-pomeron graph gets a negative AGK weight, hence 
giving in total a positive, experimentally observable cross section. However, 
the cross section for the one-pomeron cut of the two-pomeron graph is nega- 
tive. Since the one-pomeron cut of the one-pomeron graph has the same in- 
elastic final state as the one-pomeron cut of the two-pomeron graph (Fig. ^), 
one can sum both contributions to obtain a positive cross section. 

Let's consider again the example of the one- and and two-pomeron graphs. 
Concerning the topologies of the final state particles, the total cross section is 
built up of the sum of the partial cross sections of the one- and two-pomeron 
exchange graphs: (i) one-pomeron cut cti = a^^^ — 4 tT*^^-' , (ii) two-pomeron cut 
CT2 = 2 cr^^-', and (iii) diffractive cut of the two-pomeron graph CTdiff = cr^^^. 

The particle density in pseudorapidity due to a one-pomeron cut is as- 
sumed to be almost energy independent (which is true for longitudinal phase 
space models) and is denoted by dNi/drj. Provided the two-pomeron cut gives 
two times the particle yield as compared to the one-pomeron cut (in central 
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(a) 



(b) 



(c) 



Figure 9. Inelastic final states resulting from (a) the diffractive cut describing low-mass 
diflraction, (b) the one-pomeron cut, and (c) the two-pomeron cut. 



pseudorapidity region, see Fig. 
section is given by 



the inclusive inelastic charged particle cross 



da, 



ch 



dr] 



1 X (Ti 



dNi 

drj 



2 X (72 



dNi 
dr] 



X (Tdi 



iff- 



dNi 
dr] 



r(l) 



dNi 
dr] 



(23) 



Due to the topology of diffractive final states, almost no particles are produced 
in the central region in the case of a diffractive cut. Note that only the one- 
pomeron graph determines the inclusive particle cross section in the central 
region (AGK cancellation). 

On the other hand, the inclusive charged particle density follows from 



drich 



dr] 



cr(l' dNi 

CTtot dr] 



r(l) 



dNi 



r(l) 



r(2) 



dr] 



(24) 



Eq. ( P4|) allows us to understand the observed behaviour of p(0) — drich/dr] in 



pp and pp collisions. With cr^^^ ~ and ato 



„0.08 



one gets a power-law 



increase of the central particle density. This is confirmed by experimenttj: 
p(0) « 0.74 s0 i°5. 

Since the two-pomeron graph has a cross section which increases faster 
with energy than the cross section of one-pomeron graph, the model predicts 
also an increase of the multiplicity fluctuations with increasing collision enersa. 
In a model without geometric scaling this leads to violation of KNO scalingl^ 
at high energies. Furthermore, due to the characteristic structure of the one- 
and two-pomeron cuts, strong long-range correlations in pseudorapidity are 
naturally explained. For a detailed discussion, see for example □ and Refs. 
therein. 

Finally, as a representative example, model predictions for the multiplicity 
distribution and-the pseudorapidity distribution (calculated with the Phojet 
event generatorcJ) are compared with collider datall3 in Fig. |lO|. 
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UASdata i 
PHOJET 
1 cut Pomeron 

2 cut Pomerons 

3 cut Pomerons 

4 cut Pomerons 

5 cut Pomerons 




CDF 1800 GeV 
UA5 546GeV 
UA5 200GeV 
UA5 53GeV 
PHOJET 




-5 5 

pseudorapidity r| 



Figure 10. (a) Decomposition of the multiplicity distribution in pp collisions according 
to the number of generated pomeron cuts at ^/s = 200 GeV. (b) Energy-dependence of 
charged particle pseudorapidity density in pp collisions. 



5 Summary 



Regge theory, AGK cutting rules, QCD in the hmit of large numbers of col- 
ors and flavours, and the geometrical interpretation of high-energy scattering 
prove very useful in understanding the basic properties of cross sections and 
multiparticle production. They are some of the tools available today to study 
soft processes theoretically. Mailt, models combine them to obtain a detailed 
description of soft processes EIEjO'EjO. 

However we are far from being able to reliably calculate predictions for 
most soft production processes. All approaches or models discussed here have 
severe shortcomings. For example, Regge theory does not predict transverse 
momentum- related quantities, and most models implement AGK cutting rules 
without energy-momentum conservation, to name but a few. Furthermore 
the rather large number of free parameters limits the predictive power of 
calculations. 
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